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ABSTRACT:  The  interaction  of  small  molecules  (CC14,  CS2,  H20,  and  acetone)  with  single¬ 
layer  graphene  (SLG)  has  been  studied  under  steady-state  conditions  using  infrared  multiple- 
internal-reflection  spectroscopy.  Adsorption  results  in  a  broad  and  intense  absorption  band, 
spanning  the  ~200  to  500  meV  range,  which  is  attributed  to  electronic  excitation.  This  effect, 
which  has  not  previously  been  reported  for  SLG,  has  been  further  investigated  using  dispersion- 
corrected  density  functional  theory  to  model  the  adsorption  of  H20  on  SLG  supported  on  an 
Si02  substrate.  However,  the  ideal  and  defect-free  model  does  not  reproduce  the  observed 
adsorption-induced  electronic  transition.  This  and  other  observations  suggest  that  the  effect  is 
extrinsic,  possibly  the  result  of  an  adsorption-induced  change  in  the  in-plane  strain,  with 
important  differences  arising  between  species  that  form  liquid-like  layers  under  steady-state 
conditions  and  those  that  do  not.  Furthermore,  the  C— H  stretching  modes  of  CH2  groups,  incorporated  in  the  SLG  as  defects, 
undergo  nonadiabatic  coupling  to  the  electronic  transition.  This  leads  to  pronounced  antiresonance  effects  in  the  line  shapes,  which 
are  analyzed  quantitatively.  These  results  are  useful  in  understanding  environmental  effects  on  graphene  electronic  structure  and  in 
demonstrating  the  use  of  the  vibrational  spectroscopy  of  H-containing  defects  in  characterizing  SLG  structure. 


1.  INTRODUCTION 

Changes  in  the  electronic  properties  of  single-wall  carbon 
nanotubes  (SWCNTs)  and  graphene,  when  brought  into  contact 
with  molecular  vapors,  have  long  been  recognized  as  a  basis  for 
chemical  sensing1-3  and  for  doping.4'5  However,  the  mechan¬ 
isms  whereby  these  changes  occur  are  in  many  cases  not 
completely  understood.  In  many  interpretations,  the  substrate 
is  described  in  terms  appropriate  to  the  Drude  model  for  a  free- 
electron  gas.  In  this  context,  an  electron  donor  or  acceptor  (e.g., 
NH3  or  N02,  respectively)  changes  the  density  of  free  carriers 
and  thereby  the  conductivity.  For  molecules  with  a  significant 
dipole  moment,  a  decrease  in  conductivity  has  been  shown6  to 
result  from  scattering  of  free  carriers  by  the  dipolar  electric  field. 
Such  phenomena  are  well  documented7  in  the  case  of  adsorption 
on  metals.  Other  processes,  such  as  adsorbate-induced  changes 
in  the  barrier  height  at  the  SWCNT/metal  contact,  have  also 
been  implicated.8 

In  the  present  work,  it  will  be  shown  that  weak  electron 
acceptors  (CC14,  CS2,  H20)  and  donors  (acetone),  especially  in 
condensed  liquid-like  layers,  can  have  a  pronounced  effect  on  the 
electronic  structure  of  single-layer  graphene  (SLG).  This  will  be 
detected,  using  infrared  (IR)  spectroscopy  under  steady-state 
conditions,  through  observation  of  changes  in  the  broad-band 
absorption  due  to  electronic  excitation  and  of  the  coupling  of 
local  vibrational  modes  with  this  continuum.  The  line  shape  of 
the  narrow  vibrational  resonance  is  sensitive  to  damping  via 
nonadiabatic  coupling  to  the  electronic  continuum,  the  nature  of 
which  can  change  when  SLG  interacts  with  an  adsorbate.  The 
electronic  effects  reported  here  constitute  a  potential  new 
transduction  mechanism  for  chemical  sensing  as  well  as  a 
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framework  for  understanding  environmental  effects  on  the 
electronic  properties  of  SLG. 

Graphene  layers  were  deposited  on  both  sides  of  a  lightly 
oxidized  Si  parallelepiped,  configured  for  IR  multiple-internal- 
reflection  spectroscopy  (MIRS),  and  characterized  using  X-ray 
photoemission  spectroscopy  (XPS)  and  Raman  spectroscopy. 
The  MIRS  technique  permits  detection  of  surface  species  with 
high  sensitivity  under  steady-state  conditions,  which  is  necessary 
due  to  the  weak  bonding  interaction  between  SLG  and  the 
species  of  interest.  For  this  study,  the  samples  were  not  annealed, 
unless  otherwise  stated,  so  that  the  graphene  would  as  much  as 
possible  resemble  a  “free-standing”  layer.  Annealing  graphene  in 
contact  with  Si02  has  been  shown9  to  have  pronounced  effects 
on  the  chemical,  physical,  and  electronic  properties  due  to 
increased  contact  with  the  substrate. 

2.  EXPERIMENTAL  DETAILS 

2.1.  IR  Spectroscopy.  The  MIRS  data  were  obtained  using  a 
Thermo  Scientific  8700  Fourier-transform  IR  (FTIR)  system  with  a 
liquid-N2-cooled  “narrow-band”  Hg*Cdi_*Te  (MCT-A)  detector  at  a 
resolution  of  4  cm  1 .  1000  scans  were  averaged  in  about  8  min,  and 
2-fold  zero-filling  and  triangle  apodization  were  applied  to  the  inter- 
ferogram  before  transformation.  Polarization,  when  used,  was  provided 
by  a  wire-grid/KRS-5  device.  The  substrate  was  a  25  X  15  X  1  mm3 
float-zone  (FZ)  Si  parallelepiped  with  a  graphene  film  on  both  faces  and 
an  internal-reflection  angle  of  6  =  60°.  The  use  of  FZ  Si  avoids  the  strong 
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1106  cm-1  O  impurity  band  found  in  Czochralski  Si.  High-resistivity  Si 
was  used,  as  is  typical  in  IRMIRS,  to  avoid  strong  free-carrier  absorption. 
6  was  chosen  to  be  much  higher  than  the  critical  angle  (6C  ^  17°  for  Si  in 
the  mid-IR)  for  two  reasons.  First,  although  a  high  0  reduces  somewhat 
the  sensitivity  to  absorbing  thin  films,10  it  also  minimizes  the  total  optical 
path  length  through  the  Si,  thus  moving  the  transmission  cutoff  to  lower 
energy.11  Second,  a  high  6  reduces  the  possibility  of  anomalous  internal- 
reflection  phenomena10  that  can  occur  when  6  is  close  to  0C.  The 
accessible  range  was  ~  1000—7500  cm-1,  limited  on  the  low  end  by  the 
Si  transmission  and  on  the  high  end  by  the  detector  response.  Note  that 
the  spectrometer  low-pass  filter  and  the  sampling  interval  had  to  be  set 
appropriately  to  pass  Fourier  frequencies  above  the  standard  4000  cm  1 
upper  limit. 

The  substrate  was  mounted  between  two  hollow  blocks  of  Teflon  and 
the  whole  assembly  squeezed  tightly  together  to  form  a  vapor-tight  seal 
around  the  edges  of  the  Si.  There  were  1 1  internal  reflections  sampling 
graphene  not  in  contact  with  Teflon.  A  stream  of  pure  dry  N2  was  mixed 
with  another  stream  of  N2  saturated  with  the  species  of  interest  at  the 
ambient-temperature  vapor  pressure  (P0).  The  desired  partial  pressure 
(P/P0)  was  obtained  by  adjusting  the  flow  rates  of  the  two  streams, 
which  were  then  mixed  and  allowed  to  flow  continuously  over  both  faces 
of  the  sample  while  MIRS  data  were  recorded.  The  total  flow  rate  (both 
streams)  was  about  100  mL/min,  with  a  total  cell  volume  of  5.3  mL. 
Spectra  are  presented  in  the  form  of  d R/R,  the  fractional  change  in 
reflectance  per  internal  reflection  caused  by  the  introduction  of  the 
molecular  species.  The  reference  spectrum  in  all  cases  was  obtained  with 
pure  N2  flowing  through  the  cell.  No  smoothing  or  background 
correction  was  applied  to  any  of  the  data. 

Reagent-grade  CC14,  CS2,  and  acetone  and  locally  prepared  deionized 
H20  were  used  without  further  purification  and  were  contained  in  fritted 
gas-washing  bottles,  or  “bubblers”,  through  which  dry  N2  was  flowed  for 
preparation  of  the  saturated  N2  stream  mentioned  above.  To  check  for 
impurities,  the  vapor  stream  exiting  the  Teflon  cell  could  be  directed 
through  a  10  cm  gas  cell,  with  KBr  windows,  mounted  in  the  FTIR 
sample  compartment  and  the  transmission  spectrum  compared  to 
reference  data.  All  liquids  were  degassed  by  purging  with  N2.  The  liquid 
reservoirs  were  all  at  nominal  room  temperature  (25  °C)  except  for  CS2, 
which  was  maintained  at  5  °C.  The  vapor  pressures  at  these  tempera¬ 
tures  are  P0  =  113,  158,  228,  and  23  Torr,  respectively,  for  CC14,  CS2, 
acetone,  and  H20  based  on  published  parameters12  (A,  B,  and  C)  for  the 
Antoine  equation  [log10(P0)  =  A  —  B/(C  +  T)]. 

2.2.  X-ray  Photoemission  and  Raman  Spectroscopies.  The 
XPS  data  were  obtained  using  a  Thermo  Scientific  K- Alpha  system  with 
a  monochromated  Al  Ka  source  ( hv  =  1486.6  eV)  and  a  hemispherical 
electron  energy  analyzer  at  a  resolution  of  about  0.6  eV.  The  area 
sampled  was  400  /Urn  in  diameter,  and  a  flood  gun  was  used  for  sample 
neutralization  because  the  Si  substrates  used  for  MIRS  are  highly 
resistive.  The  BE  was  determined  using  a  reference  value13  of  99.34 
eV  for  the  elemental  Si  2p3/2  peak.  The  BE  correction  required  was 
small,  typically  ~0.1  eV.  The  C  Is  least-squares  fits  were  done  using  a 
Gaussian -broadened  Doniach— Sunjic  line  shape  for  the  graphene  peak 
and  pure  Gaussians  for  other  features,  if  any.  These  were  added  to  a 
polynomial  background  function,  and  all  parameters  pertaining  to  the 
peak  and  the  background  were  unconstrained  in  the  fitting  process. 

Raman  spectra  were  acquired  using  confocal  Raman  system  with  a 
single-mode  532  nm  laser  and  a  Princeton  Instruments  CCD  array. 

2.3.  Sample  Preparation.  Graphene  films  were  grown  via  low- 
pressure  chemical  vapor  deposition  (CVD)  on  Cu  foils.14  Twenty-five 
micrometer  thick  Cu  foils  (Alfa  Aesar,  99.8%)  were  heated  to  1030  °C 
under  flowing  H2  at  ~600  mTorr  (with  a  system  base  pressure  of  ~5 
mTorr).  At  the  growth  temperature,  the  H2  pressure  was  decreased  to 
<20  mTorr,  and  CH4  was  introduced  at  pressures  ranging  from  40  to 
300  mTorr  and  flow  rates  from  5  to  80  seem  for  15—45  min,  after  which 
the  sample  was  quenched  to  room  temperature.  Methane  pressures  and 


growth  times  were  varied  to  adjust  the  density  of  bilayer  islands 
from  <~4%  to  approximately  40%  coverage. 

In  another  experiment,  we  attempted  to  introduce  deuterium  (D) 
into  the  films  by  flowing  high-purity  D2/ CD4  during  growth  under  the 
same  conditions  described  above  for  H2/ CH4.  However,  as  described  in 
the  main  text,  these  films  did  not  show  evidence  of  D  in  the  IR  MIRS 
experiments.  We  speculate  this  may  be  due  to  (i)  the  presence  of 
nondeuterated  hydrocarbons  during  growth,  (ii)  D  exchanging  with  H 
during  the  etching  and  transfer  process  (described  below),  or  (iii) 
possibly  atmospheric  H20  decomposing  at  defect  sites  in  the  film.  The 
most  likely  location  for  C— H  bonds  in  CVD  graphene  is  at  domain 
boundaries,  vacancies,  or  free  edges  (such  as  those  associated  with 
isolated  multilayer  spots). 

After  growth,  one  side  of  the  graphene-covered  Cu  foil  was  coated 
with  poly(methyl  methacrylate)  (PMMA),  and  then  the  Cu/graphene/ 
PMMA  sample  was  floated  in  a  Transene  Cu  etchant.  Afterward,  the 
PMMA/ graphene  film  was  transferred  to  a  water  bath.  Films  prepared  in 
this  manner  were  subsequently  transferred  onto  both  sides  of  a  Si  MIRS 
parallelepiped  and  then  soaked  in  acetone  to  remove  the  PMMA  For 
two  of  the  samples,  post-transfer  annealing  was  carried  out  in  a  25  mm- 
diameter  tube  furnace  with  flowing  Ar  at  atmospheric  pressure  at  either 
T  =  150  °C  or  T  =  400  °C  for  1  h. 

3.  RESULTS  AND  DISCUSSION 

3.1.  Sample  Characterization.  Figure  la  shows  C  Is  XPS 
data,  which  can  be  least-squares  fitted  with  a  single  Gaussian- 
broadened  Doniach— Sunjic  (DS)  function.15  The  DS  theory 
describes  the  screening  of  a  photoionized  core  hole  in  a  free-electron 


Figure  1.  C  Is  XPS  data  for  a  graphene  sample  (a)  with  no  detectable  O 
contamination  (sample  Gl)  and  (b)  with  a  small  concentration  of  sp3 
C— O  bonds  (sample  G2).  The  points  are  the  raw  data  after  removal  of 
the  fitted  polynomial  background  function.  The  lines  are  least-squares 
fits  of  a  Gaussian-broadened  DS  function  for  the  graphene  peak  and  (in 
panel  b)  a  pure  Gaussian  for  the  C— O  satellite.  The  residuals  (fit  minus 
data,  not  statistically  weighted)  are  shown  multiplied  by  a  factor  of  10. 
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metal,  which  results  in  a  distorted  Lorentzian  line  shape  with  a 
tail  extending  to  higher  binding  energy  (BE).  The  DS  function 
is  only  approximately  accurate  in  the  case  of  graphite,16  which 
is  a  semimetal.  Nevertheless,  the  fit  shows  only  sp2-hybridized 
C  (BE  =  284.5  eV)  with  little  or  no  evidence  of  sp3-hybridized 
C  or  of  C  — O  bonds  (as  in  C  — O  — C  or  C  — OH  groups).  Such 
species  would  result  in  satellites17  higher  in  BE  by  about  0.5  or 
1.8  eV,  respectively,  than  the  sp2  C  Is  peak.  The  DS  asymmetry 
parameter  (a  =  0.13)  agrees  very  well  with  the  result16  (a  = 
0.125)  for  bulk  highly  oriented  pyrolytic  graphite  (HOPG). 
This  sample,  henceforth  referred  to  as  “Gl”  (for  “graphene 
sample  no.  1”),  exhibited  a  C  Is  intensity  that  was  spatially 
uniform  (±10%)  over  most  of  both  sides  of  the  MIRS 
substrate,  indicating  a  minimal  density  of  holes  (i.e.,  exposed 
Si02). 

Another  sample  (G2,  Figure  lb)  showed  a  weak  satellite  at 
1.90  eV  higher  BE,  with  an  integrated  intensity  equal  to  a  few 
percent  of  that  of  the  graphene  peak,  which  indicates  a  small 
concentration  of  sp3  C— O  bonds.  There  is  also  evidence  in  the 
residual  spectrum  of  structure  near  288  eV  due  to  C=0  bonds.17 
Some  areas  showed  about  the  same  C  Is  intensity  as  did  sample 
Gl,  but  others  showed  a  much  smaller  intensity,  indicating  holes 
in  the  graphene  layer.  A  slightly  larger  DS  asymmetry  (a  =  0.17) 
was  needed  to  fit  the  graphene  peak,  which  suggests  the  presence 
of  a  weak  and  unresolved  sp3  component  at  slightly  higher  BE. 
Comparison  of  data  for  samples  with  different  degrees  of 
perfection  will  be  useful  in  the  following  discussion.  In  all,  five 
different  SLG  and  three  multilayer  graphene  samples  were 
studied,  and  two  of  the  SLG  samples  were  subjected  to  annealing 
(to  be  discussed  later),  giving  in  effect  a  total  of  10  samples.  The 
XPS  data  in  Figure  1  are  representative  of  all  of  these  samples. 
Analysis  of  the  Si  2p  XPS,  described  in  section  1  of  the 
Supporting  Information,  indicates  an  oxide  layer  underlying 
the  graphene  with  a  thickness  in  the  range  of  19—28  A  for  all 
samples  studied  here.  This  layer  is  composed  largely  of  Si02  with 
little  or  no  indication  of  suboxide  features.  The  C  Is  XPS  data  are 
discussed  further  in  the  Supporting  Information,  section  2.  Two 
of  the  samples  also  showed  a  trace  amount  of  Cu  remaining  from 
growth,  but  the  Cu/C  atom  ratio  was  <0.4%.  Additional  sample 
characterization  using  Raman  spectroscopy  and  optical  micro¬ 
scopy  is  discussed  in  the  Supporting  Information,  section  3. 

3.2.  IR  Spectroscopy  Results  —  Overview  and  Electronic 
Excitation.  In  this  section,  the  general  IR-spectroscopic  results 
are  described,  and  the  adsorption-induced  electronic  excitation  is 
examined.  Figure  2  shows  MIRS  data  (sample  Gl)  for  acetone, 
H20,  and  CC14  (all  at  P/P0  =  1.0,  where  P0  is  the  room- 
temperature  vapor  pressure)  over  a  broad  range  from  1400 
to  7000  cm-1  (170—870  meV),  which  includes  the  strong 
l'(C=0)  carbonyl  stretching  mode  of  acetone  (at  1715  cm*1) 
and  the  v(OH)  stretching  and  (5(HOH)  bending  modes  of  H20 
(at  ~3300  and  1630  cm  ,  respectively).  CC14  has  no  funda¬ 
mentals  in  this  range.  Note  that  dR/R  is  defined  such  that 
downward-pointing  peaks  correspond  to  increased  absorption 
due  to  the  presence  of  the  molecular  layer.  The  v(C=0)  peak 
appears  at  1715  cm*1,  which  is  characteristic  of  liquid  acetone 
(vs  1738  cm  1  for  the  vapor),  indicating  a  condensed  multilayer. 
The  v(OH)  region  ofH20  shows  evidence  of  overlapping  bands 
at  about  3380  and  3220  cm  ,  the  latter  suggesting  an  ice-like 
ordering18'19  at  the  graphene/H20  interface.  This  aspect  is 
beyond  the  scope  of  the  present  work,  but  the  v(OH)  structure 
is  shown  and  discussed  in  more  detail  in  the  Supporting 
Information,  section  4. 


Figure  2.  MIRS  data  for  graphene  (sample  Gl)  exposed  to  acetone, 
H20,  and  CC14  (all  at  P/P0  =  1.0).  All  spectra  are  shown  on  the  same 
dR/R  scale  (normalized  to  the  total  number  of  internal  reflections)  and 
without  any  arbitrary  vertical  displacement.  The  dashed  line  shows  the 
position  of  dR/R  =  0.  The  increased  noise  at  the  high-energy  end  is  due 
to  decreasing  detector  sensitivity.  The  intense  acetone  l'(C— O)  peak  at 
1715  cm  1  has  been  truncated  for  clarity  of  display.  The  structure  at  the 
low-energy  end,  near  1450  cm  1 ,  results  from  a  strong  multiphonon 
absorption  peak  in  the  Si  substrate  that  does  not  completely  cancel  in 
dR/R.  The  weak  feature  labeled  “l'(C— O)”  is  due  to  an  impurity 
discussed  in  the  text.  Note  that  an  increase  in  absorption  corresponds  to 
downward-pointing  peaks. 


The  weak  feature  at  the  extreme  low-energy  end  of  all  spectra 
("-'1450  cm  *)  is  due  to  a  strong  Si  multiphonon  absorption  that 
does  not  completely  cancel  in  the  ratio  of  single-beam  spectra.  A 
very  weak  structure,  most  noticeable  in  the  p-polarized  CC14 
spectrum,  is  seen  near  1700  cm*1  and  is  thought  to  arise  from  a 
perturbation  of  the  v(C=0)  mode20  of  carbonyl  impurities  by 
interaction  with  the  liquid  multilayer.  Although  such  species  are 
at  or  below  the  XPS  detection  limit  (~1%  of  the  graphene  C  Is 
intensity,  cf.,  Figure  l),  they  may  be  detectable  in  MIRS  due  to 
the  large  oscillator  strength  of  this  mode.  However,  there  is  no 
indication  of  the  v(C=C)  stretching  mode,  which,  for  defect- 
free  SLG,  is  IR-forbidden21  and  not  detected  in  transmission 
spectra.20  This  mode,  which  produces  the  Raman-active  G-band 
found  at  1588  cm  1  in  SLG,  will  be  discussed  later  in  connection 
with  MIRS  data  for  few-layer  graphene  (FLG). 

The  features  of  interest  in  the  present  work  are  the  very  intense 
and  broad  absorption  bands  spanning  the  ~H600— 4000  cm*1 
(200—500  meV)  range  and  the  unusual  line  shapes  seen  in  the 
v(C— H)  stretching  region  (~2850— 2950  cm*1).  Purging  with 
dry  N2  at  room  temperature,  following  exposure  to  any  of  the 
species  studied  here,  eventually  removes  these  features.  This 
indicates  that  they  are  not  due  to  covalent  bonding  between  the 
reagents  and  either  graphene  or  the  Si02  underlayer.  Rather,  these 
effects  are  reversible  at  room  temperature  and  are  observable  only 
under  steady-state  conditions.  Numerical  simulation  of  the  MIRS 
experiment  (Supporting  Information,  section  5)  shows  that  they 
do  not  arise  from  thin-film  optical  effects,22  and  none  of  these 
effects  were  observed  in  “blank”  experiments  for  which  no 
graphene  was  deposited  on  the  oxidized  Si  surface  (Supporting 
Information,  section  6).  The  blank  sample  was  processed  using  the 
same  procedure  as  for  the  graphene  samples  including  the  use,  and 
subsequent  removal,  of  PMMA  as  described  above  (section  2.3). 
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Further  observations  regarding  the  blank  sample  are  given  in  the 
Supporting  Information,  sections  2  and  6. 

Data  obtained  in  the  range  of  0.1  <  P/P0  <  1.0  (not  shown) 
indicate  absorption  intensities  that  increase  with  P/P0.  The 
absorption  spectrum  of  ideal  SLG  does  not  exhibit  any  structure 
in  this  energy  region.23  It  has  been  shown  that  movement  of  the 
Fermi  level  (£p)  in  bilayer  graphene  (BLG)  or  in  FLG  by 
electron  transfer  to  or  from  an  adsorbate24  or  by  application  of 
an  electric  field25-28  can  change  the  nature  of  the  interband 
spectrum  in  the  mid-IR.  However,  the  interband  transitions  in 
BLG  and  FLG,  which  are  modified  by  these  perturbations,  are 
not  found  in  pristine  SLG.  Hence,  the  mechanisms  discussed  in 
refs  24—28  are  not  believed  to  be  directly  applicable  to  the 
present  data.  A  subsequent  section  will  discuss  the  effect  of 
doping  on  the  MIRS  data;  however,  doping  (i.e.,  a  shift  in  £P) 
cannot  by  itself  account  for  the  appearance  of  a  mid-IR  absorp¬ 
tion  band  in  SLG.  The  Supporting  Information,  section  7, 
discusses  similar  data  for  sample  G2,  which  has  known  imperfec¬ 
tions  as  mentioned  above.  The  spectra  are  quantitatively  differ¬ 
ent,  although  broad-band  absorptions  and  v(C— H)  line  shape 
anomalies  are  again  seen  for  all  three  adsorbates.  For  all  samples, 
the  intensity  of  the  broad-band  absorption  consistently  varies  as 
CC14  <  H20  <  acetone. 

Numerical  evaluation  of  the  IR  electric  field  throughout  the 
MIRS  sample  (Supporting  Information,  section  8)  shows  that 
|£pjZ|2,  the  intensity  of  the  surface-normal  component  of  the 
p-polarized  field,  is  small  within  the  graphene  layer  and  that  the 
in-plane  field  intensities,  |£pi7|2  and  |£sa.|2,  are  relatively  large  and 
nearly  equal.  Hence,  the  dynamic  dipole  moment  associated  with 
the  electronic  excitation  must  lie  in  the  SLG  plane.  Given  the 
symmetry  of  ideal  SLG,  one  then  expects  the  s-  and  p-polarized 
spectra  to  be  identical,  which  is  not  the  case.  This  could  be  an 
indication  of  in-plane  strain,  which  is  found  in  theoretical 
work29,30  to  have  a  pronounced  effect  on  the  electronic  structure 
of  SLG.  For  certain  directions,  the  effect  of  strain  on  the  density 
of  states  (DOS)  is  particularly  apparent  near  £p,  within  the 
energy  range  studied  here.  Furthermore,  the  extent  of  the  IR 
polarization  dependence  seen  experimentally  varies  from  sample 
to  sample,  which  is  consistent  with  an  extrinsic  effect  such  as 
strain.  Hence,  a  change  in  strain  due  to  contact  with  the  liquid- 
like  layer  is  tentatively  suggested  as  a  possible  mechanism  giving 
rise  to  the  electronic  absorption  band.  Buckling,  or  corrugation,  is 
another  well-known  extrinsic  effect  in  SLG  supported  on  Si02 
(e.g.,  refs  9,31).  However,  this  is  found32  not  to  affect  the  DOS 
near  £F  except  in  extreme  cases. 

3.3.  Computational  Model.  In  this  section,  we  consider 
whether  the  observed  electronic  absorption  can  be  explained  in 
terms  of  ideal  SLG  supported  on  Si02,  using  H20  as  the  simplest 
adsorbate  and  one  which  is  observed  (Figure  2)  to  give  a 
pronounced  effect.  There  is  also  a  body  of  prior  theoretical  work 
for  H20/SLG  with  which  to  compare  the  present  results.  It  has 
previously  been  proposed33,34  that  dangling-bond  defects  on  the 
Si02  surface  are  important  in  the  response  to  adsorbates  for  SLG 
supported  on  Si02.  However,  any  such  highly  reactive  defects  on 
air-exposed  Si02  would  be  passivated  by  dissociative  adsorption 
of  H20  vapor  to  form  OH  groups  and  would  not,  therefore,  be 
expected  at  a  real  graphene— Si02  interface. 

Some  density  functional  theory  (DFT)  treatments  of  H20  on 
free-standing  SLG  have  employed  methods  that  do  not  com¬ 
pletely  account  for  dispersion  (i.e.,  noncovalent  or  van  der  Waals 
interaction),  which  is  expected  to  be  a  major  factor  in  this  system. 
The  adsorption  energies  (A£ads)  found  in  these  studies,33"’5 


(a) 


Figure  3.  Computational  results  for  the  structure  of  H20  on  (a)  free¬ 
standing  SLG  and  on  (b)  SLG  on  Si02.  In  (b),  a  7-layer,  (lll)-oriented 
slab  is  cut  from  the  bulk  /3-cristobalite  lattice  and  the  under-coordinated 
Si  atoms  terminated  with  OH.  A  (2  X  2)  SLG  supercell  is  placed  on  top, 
and  one  H20  per  supercell  is  added  on  top  of  the  SLG.  The  atom  colors 
are  off-white  (H),  black  (C),  red  (O),  and  beige  (Si).  See  text  and  the 
Supporting  Information,  section  9,  for  further  details. 


using  “pure”  generalized  gradient  approximation  functionals, 
have  been  <40  meV/H20  with  little  or  no  effect  on  the  graphene 
electronic  structure  near  £F  for  coverages  up  to  a  monolayer.  On 
the  other  hand,  studies  using  coupled-cluster  methods36  or 
symmetry-adapted  perturbation  theory,37  which  accurately  treat 
dispersion,38  find  A£ads  of  130—135  meV/H20.  Other  work,’9 
using  dispersion-corrected  DFT,  finds  a  smaller  A£ads  but  one 
that  is  still  larger  than  when  noncovalent  interaction  is  neglected. 
These  studies  all  agree  that  in  the  lowest-energy  configuration 
(Figure  3a)  H20  is  oriented  normal  to  the  surface  with  the  H 
atoms  downward  (i.e.,  closer  than  O  to  the  graphene).  However, 
the  effects,  if  any,  on  the  SLG  electronic  structure  were  not 
reported. 

We  have  performed  DFT  calculations  similar  to  those  de¬ 
scribed  in  ref  33  but  employing  methods  that  include  dispersion. 
Full  details  are  given  in  the  Supporting  Information,  section  9. 
The  goal  is  to  determine  whether  dispersion,  which  increases 
A£ads,  also  affects  the  SLG  electronic  structure  near  £p.  A 
coverage  of  0.25  H20  per  SLG  unit  cell  was  used,  which 
(Supporting  Information,  section  9)  is  close  to  the  coverage  of 
“terminal”  H20  on  the  (0001)  surface  of  hexagonal  ice.  Thus,  the 
calculation  simulates,  in  effect,  SLG  in  contact  with  a  hypothe¬ 
tical  monolayer  of  ice.  For  free-standing  SLG,  a  A£ads  of  207 
meV/H20,  not  corrected  for  basis  set  superposition  error 
(BSSE),  was  found.  The  lack  of  a  BSSE  correction,  which  would 
reduce  the  magnitude  of  AEads,  makes  it  difficult  to  compare  this 
result  with  previous  reports. 36,37,39  However,  without  dispersion, 
no  binding  (i.e.,  A£ads  0)  was  found,  confirming  the  impor¬ 
tance  of  including  this  effect.  In  this  preliminary  study,  we  have 
not  determined  the  effect  of  higher  H20  coverages.  Even  with 
dispersion,  no  effect  due  to  H20  was  seen  in  the  DOS  (not 
shown)  within  ~1.5  eV  of  £P.  Also,  no  doping  was  observed,  in 
the  form  of  a  shift  of  £p  into  the  valence  or  conduction  band, 
because  this  appears33  only  for  an  H20  multilayer  with  ice-like 
ordering. 
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Figure  4.  Similar  to  Figure  2  but  showing  the  r(C— H)  region.  The  dR/R 
scale  applies  to  the  acetone  and  CC14  data.  The  H20  spectra  have  been 
multiplied  by  a  factor  of  3,  and  different  spectra  have  been  displaced 
vertically  for  clarity. 


Following  ref  33,  Si02  was  included  using  a  slab  (Figure  3b) 
derived  from  the  (ill)  surface  of  /3-cristobalite  with  dangling  Si 
bonds  terminated  in  OH.  The  lattice  constant  of  /5-cristobalite  is 
fl0  -  7.16  A;  hence,  that  of  the  (ill)  surface  unit  cell  (a0/  2)  is 
only  '•“2.7%  larger  than  twice  that  of  graphene.  Thus,  a  (2  x  2)  SLG 
supercell  is  commensurate  with  a  slightly  compressed  /3-cristobalite 
(ill)  surface.  Previous  work33  (omitting  dispersion)  for  SLG 
interacting  with  defect-free  Si02  reported  no  significant  effects  in 
the  vicinity  of  the  SLG  Fermi  level,  and  the  present  results  concur 
with  this  finding.  The  present  work  allows  an  estimate  of  the  energy 
of  adhesion  between  SLG  and  Si02  as  E adh  =  E(Si02)  +  £(SLG)  — 
E( Si02  +  SLG),  where  E  is  the  total  energy  after  geometry 
relaxation.  The  result  is  £adh  =  320  meV  (not  corrected  for  BSSE), 
which,  upon  normalization  to  the  area  of  the  (2  X  2)  SLG 
supercell,  gives  1.5  eV/nm2.  This  can  be  compared  to  a  value  of 
about  0.50  eV/nm2  (80  mj/m2)  obtained40  for  vitreous  Si02  and 
graphite  using  molecular  mechanics  with  a  parametrized  force  field. 

The  next  step  is  to  add  H20  to  the  Si02-supported  SLG.  In  the 
relaxed  structure  (Figure  3b),  H20  lies  nearly  flat,  with  AEads  = 
257  meV/H20.  The  difference  of  50  meV/H20  in  A Eads  from 
that  for  free-standing  SLG  given  above  should  be  essentially 
unaffected  by  BSSE.  Thus,  the  underlying  Si02  slightly  enhances 
the  adsorption  of  H20;  however,  there  is  again  no  effect  on  the 
DOS  near  Ep.  It  is  concluded  that  pristine  SLG  is  weakly 
hydrophilic  but  that  the  MIRS  results  cannot  be  explained  simply 
in  terms  of  an  ideal  H20/SLG/Si02  structure  and  that  some 
extrinsic  effect,  such  as  strain,  is  needed  to  account  for  the  data. 

3.4.  IR  Spectroscopy  Results  —  Vibrational  Line  Shapes 
and  Fine  Structure.  We  now  proceed  to  the  anomalous  v(C— H) 
line  shapes.  Figure  4  shows  the  same  data  as  in  Figure  2  but 
focusing  on  the  C— H  stretching  region.  All  three  adsorbates  gave 
similar  features  near  2920  and  2850  cm-1,  although  they  are 
relatively  weak  in  the  case  of  H20.  A  reason  for  the  small  effect  of 
H20  will  be  suggested  in  the  following  section.  These  features  are 
reproducible  for  all  10  samples,  and  spectra  obtained  at  2  cm-1  (vs 
4  cm-1)  resolution  showed  a  negligible  reduction  in  line  width, 
indicating  that  the  data  are  not  resolution  limited.  The  higher- 
energy  feature  is  essentially  independent  of  polarization,  whereas 
the  lower-energy  component  is  stronger  in  p-polarization. 


-i 

cm 


Figure  5.  (Top)  Data  obtained  under  steady-state  exposure  to  (a)  CC14 
and  (b)  CS2  at  P0  =  99  Torr  (22  °C)  and  158  Torr  (5  °C),  respectively, 
with  an  unpolarized  IR  beam.  The  data  were  obtained  for  sample  G2, 
and  the  spectra  have  been  displaced  vertically  for  clarity.  The  weak 
structure  near  3700  cm  1  is  a  baseline  artifact.  (Bottom)  Least-squares 
fit  of  eq  1  to  an  s-polarized  CC14  spectrum.  The  points  are  the  raw  data, 
after  subtraction  of  the  fitted  polynomial  background,  and  the  solid  line 
is  the  fit.  Note  the  different  energy  scales  in  the  upper  and  lower  panels. 
The  fitting  parameters  are:  v0  =  2918  cm  1 ,  v0t  =  —3.49,  TG  =  4.6  cm  ', 
TL  =  4.7  cm  1'l  and  v0  =  2844  cm  *,  v0t  =  —  0.01,  TG  =  5.3  cm  ',  T|  = 
4.2  cm  ' .  The  inset  shows  a  configuration  of  the  CH2  that  is  consistent 
with  the  different  v0z  results  (see  text).  The  dashed  lines  show  the  IR 
beam,  and  the  s-  and  p-vectors  define  the  polarizations.  and 

depict  the  vibrational  and  electronic  DDMs.  The  rotation  of  the  CH2 
about  the  surface  normal  is  arbitrary;  hence,  vs  can  be  excited  only  by  the 
surface-normal  component  in  p-polarization,  but  l'as  can  be  excited  by 
s-polarization  and  by  the  in-plane  component  in  p-polarization. 


However,  in  the  case  of  H20,  the  very  weak  v(C— H)  structure 
does  show  a  polarization  dependence,  being  almost  undetectable 
in  s-polarization.  Additional  structure  is  seen  for  acetone  near 
3000  cm-1  due  to  the  molecular  v(C— H)  modes.  It  should  be 
noted  that  the  graphene  was  not  intentionally  hydrogenated; 
hence,  the  v(C— H)  modes  are  due  to  defects  and/or  edge  sites. 

The  following  discussion  will  focus  on  the  CC14  results  because 
CCI4  presents  no  interfering  absorption  bands  in  the  region  of 
interest  and  is  “inert”  in  the  sense  of  having  no  dipole  moment 
and  only  a  weak  interaction  with  graphene-like  materials.41 
Figure  5  compares  data  obtained  during  exposure  to  CC14  and 
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to  CS2;  another  small,  nonpolar  molecule  with  no  C— H  bonds. 
The  results  are  clearly  very  similar.  These  data  were  obtained  for 
sample  G2  and  are  in  good  agreement  with  the  CC14  data  for 
sample  G1  (Figure  2).  As  mentioned  above  (section  2.3), 
unsuccessful  attempts  were  made  to  incorporate  D  using  CD4 
and  D2  in  place  of  CH4  and  H2  in  the  growth.  However,  it  was 
later  determined  that  nondeuterated  hydrocarbons  were  also 
present  in  the  growth  ambient  and  were  possibly  the  source  of 
the  persistent  CHX  species.  Hence,  it  cannot  yet  be  determined  at 
what  point  in  the  sample  growth  and  preparation  the  C— H 
bonds  are  formed.  Because  the  CHX  groups  involve  sp3  sites,  a 
firm  upper  limit  of  ~0.1  is  estimated  for  the  CH^/total-C  atomic 
ratio,  based  on  the  maximum  possible  sp3/sp2  ratio  in  the  C  Is 
XPS  (Supporting  Information,  section  2). 

The  energies  of  the  2850  and  2920  cm-1  peaks  are  consistent 
with,  respectively,  the  symmetric  (in-phase,  l's)  and  asymmetric 
(out-of-phase,  vas)  C— H  stretching  modes42  of  an  sp3-hybri- 
dized  CH2  group,  which  would  be  a  defect  in  SLG.  The  line  shape 
of  the  vas(CH2)  mode  is  unusual,  being  distorted  and  inverted  or 
derivative-like  in  appearance.  This  effect  is  much  less  pro¬ 
nounced,  and  may  even  be  absent,  in  the  vs(CH2)  mode.  The 
anomalous  line  shape  results  from  the  coupling  of  the  vibrational 
mode  to  the  electronic  excitation  continuum.  The  theory  for  this 
effect  (termed  “antiresonance”,  “vibronic  interference”,  or  “the 
Fano  effect”)  was  developed  in  the  context  of  adsorbate  vibra¬ 
tional  spectra  by  Langreth  et  al 43,44  and  reviewed  from  an 
experimental  perspective  by  Chabal.45  The  dynamic  dipole 
moment  (DDM,  fl)  is  the  sum  of  a  nuclear  (/tn)  and  an  electronic 
(jle)  component.  In  the  case  of  coupling,  fln  and  fle  are  out  of 
phase  resulting  in  a  complex  DDM  given  by  fi  =  fi1  +  ifi"  =  //  ( 1  + 
tor),  which  leads  to  an  asymmetric  line  shape. 

Antiresonance  in  the  IR  spectra  of  adsorbates  has  been  seen 
for  H  on  W(l00)  and  Mo(lOO)  where  the  first  overtone  of  the 
metal-H  bending  mode  couples  to  a  broad  continuum  of  surface- 
state  excitations.46  It  has  also  been  seen  for  CO  on  Cu(lOO), 
where  the  frustrated  rotation  of  the  CO  couples  to  the 
continuum  of  adsorbate-induced  surface-state  excitations.47 
Antiresonance  has  been  observed”4  for  the  1589  cm  1  in-plane 
(Eiu)  v(C — C)  phonon  in  FLG,  where  the  strength  of  the 
electron— phonon  coupling  was  varied  using  adsorption  of  an 
electron  donor  (NH3)  to  shift  £P.  The  shift  in  £P  induces 
interband  transitions  that  can  then  couple  to  the  Elu  phonon. 
A  similar  effect  has  been  observed25-27  by  applying  an  electric 
field  to  BLG  to  open  a  band  gap,  in  which  case  the  strength  of  the 
coupling  varies  as  the  gap  is  increased  to  coincide  with  the 
phonon  energy.  The  present  result  is,  however,  fundamentally 
different  from  those  of  previous  studies  because  it  involves  SLG 
interacting  with  molecules  having  relatively  little  (vs  NH3  or 
N02)  electron  donor  or  acceptor  character.  Another  form  of 
coupling  between  electronic  and  vibrational  excitations  has  been 
seen20  in  reduced  graphene  oxide  (RGO).  This  is  different 
from  the  effect  studied  here  because  in  the  case  of  RGO  there 
is  no  phase  difference  between  the  vibrational  and  electronic 
excitations  and,  therefore,  no  antiresonance.  In  RGO,  the 
electronic  absorption  results  from  the  displacement  of  the  O 
atoms  in  the  ^(C— O)  mode  and  is  therefore  in  phase  with  the 
atomic  motion. 

In  all  of  the  examples  noted  above,  as  well  as  in  the  present 
study,  antiresonance  involves  an  in-plane-polarized  vibrational 
mode  coupling  to  an  in-plane-polarized  electronic  transition. 
Coupling  is  seen  here  (and  in  refs  20  and  24—27)  only  for 
vibrational  modes  involving  species  within  the  graphene  layer 


(e.g.,  C=C  bonds  or  impurity  C— O  or  C— H  sites)  and  not  for 
those  involving  adsorbates.  Thus,  the  C  atom  involved  in  the 
vibrational  mode  must  be  a  graphene  site  for  strong  coupling  to 
occur.  For  example,  Liu  et  al.24  observe  a  three-peaked  v(C— H) 
structure  in  the  2850—2950  cm-1  range,  which  is  characteristic 
of  an  aliphatic  hydrocarbon42  (possibly  adsorbed  from  the 
ambient)  and  which  does  not  exhibit  antiresonance.  Further¬ 
more,  none  of  the  acetone  molecular  modes  seen  in  the  present 
work  exhibit  any  line  shape  anomalies.  Thus,  the  internal 
vibrational  modes  of  the  weakly  bound  molecular  adsorbates 
studied  to  date  do  not  themselves  couple  strongly  to  graphene 
electronic  transitions.  It  is  possible,  however,  that  a  strongly 
chemisorbed  species  would  exhibit  antiresonance  as  in  the 
examples  of  H  and  CO  on  metals  cited  above. 

Figure  5  shows  a  fit  of  representative  s-polarized  CC14  data 
with  the  antiresonance  line  shape  function  for  an  isolated 
oscillator  derived  by  Crljen  and  Langreth,44  which  (with  ft)  = 
2  Jtv)  is  given  by45 

I'o/t'2  vt2  vTl(  1  -  v2r2)  -  2vr(v2  -  v2) 

Imam  = -  — — | - — -= - r - 

n  rL  ( v 2  —  Vg)  +  v2rl 

(i) 

and  which,  for  T  =  0,  reduces  to  the  Lorentzian  line  shape.  Here, 
a(v)  is  the  local  polarizability,  Tl  is  the  Lorentzian  line  width, 
and  all  other  terms  have  been  defined  previously.  An  s-polarized 
spectrum  was  used  because  the  relative  intensity  of  the  lower- 
energy  vs  peak  is  weaker  than  in  p-polarization,  which  simplifies 
the  fit.  The  fitting  is  made  difficult  by  the  presence  of  overlapping 
weak  and  poorly  resolved  structure  due  to  overtones  of  the 
(5(C— H)  bending  modes42  and  to  the  stretching  modes  of  other 
C— H  sites,  which  are  not  included  in  the  fit.  Evidence  for  the 
existence  of  other  C— H  sites  will  be  given  below.  Cl(v )  was 
convoluted  with  a  Gaussian  (full-width  at  half-maximum, 
FWHM,  =Tg)  to  include  inhomogeneous  broadening  and  finite 
spectrometer  resolution,  both  approximated  as  Gaussians.  To 
this  was  added  a  cubic  or  quartic  polynomial  background 
function,  and  all  terms  were  unconstrained  in  the  fit. 

The  Tg  values  (given  in  the  figure  caption)  are  about  the  same 
for  either  peak,  as  are  the  two  TL  values.  Subtracting  the  FTIR 
resolution  (4  cm  :)  in  quadrature  from  TG  gives  an  inhomoge¬ 
neous  broadening  of  about  2.3  (3.5)  cm  1  for  the  vas  (vs)  mode. 
The  main  difference  is  in  v0T,  a  dimensionless  asymmetry 
parameter  that  reflects  the  strength  of  the  coupling  between 
the  electronic  and  vibrational  excitations  and  that  is  large  for  the 
2918  cm  1  vas  mode  and  small  for  the  2844  cm  1  vs.  Fits  to 
other  data  sets  give  somewhat  different  results  for  v0T,  but  the 
value  for  vas  (vs)  is  always  in  the  range  of  —2.0  to  —3.5  (0  to 
—  0.5)  .  For  comparison,  v0T  in  the  range  of  —1.6  to  +0.54  has 
been  found46,47  for  adsorbates  on  metals.  The  Supporting 
Information,  section  10,  shows  how  the  line  shape  changes  over 
a  range  of  v0T  values  and  how  these  shapes  replicate  those  seen  in 
Figure  4. 

The  DDM  of  vs  is  parallel  to  the  2-fold  (C2)  axis  of  the  CH2 
group,  whereas  that  of  ras  is  perpendicular  to  the  C2  axis  and  lies 
in  the  CH2  plane.  Figure  5(inset)  shows  a  configuration  in  which 
the  CH2  bridges  two  C  atoms  and  lies  above  the  SLG  plane  and 
which  is  consistent  with  the  fitting  results.  It  is  not  clear  whether 
the  CH2  exists  as  a  strained,  cyclopropane-like  structure  or  forms 
as  a  result  of  buckling  of  the  graphene  plane  in  response  to  the 
breaking  of  C=C  bonds.  Here,  the  vs  DDM  is  normal  to  the 


11031 


dx.doi.org/10.1021/la201669j  | Langmuir  2011, 27, 1 1026-1 1036 


Langmuir 


ARTICLE 


surface  and  therefore  unable  to  couple  to  fle,  whereas  the  vas 
DDM  can  interact  strongly.  Deviations  from  this  geometry  in 
which  the  vs  DDM  has  a  finite  projection  on  the  SLG  plane  will 
increase  (decrease)  the  magnitude  of  v0T  for  vs  (vas).  In  this 
model,  vs(CH2)  is  completely  p-polarized.  However,  the  ob¬ 
served  vJVas  relative  intensity  is  about  1/3  (2/3)  in  s-  (p-) 
polarization,  not  corrected  for  the  relative  magnitudes  of  the 
different  components  of  the  evanescent  electric  field  (Supporting 
Information,  section  8).  The  fact  that  vs(CH2)  is  not  completely 
p-polarized  (i.e.,  vjv^  is  nonzero  in  s-polarization)  constitutes 
an  apparent  violation  of  the  selection  rule  for  this  geometry.  This 
could  result  from  the  buckling7’1  (or  corrugation)  of  SLG 
supported  on  Si02,  which  would  tip  the  C2  axis  away  from  the 
z-axis,  as  described  in  the  Supporting  Information,  section  11. 

A  recent  study,48  using  sum-frequency  generation  (SFG), 
reported  the  vibrational  spectrum  of  individual  H  atoms  chemi¬ 
sorbed  as  sp3  C— H  groups  on  pristine  SLG  grown  in  situ  in  an 
ultrahigh  vacuum  on  an  Ir(lll)  surface.  Because  the  s-polarized 
IR  electric  field  is  suppressed  at  the  vacuum/metal  surface,45  only 
modes  with  a  DDM  normal  to  the  surface  were  detected.  These 
were  superimposed  on  a  broad  continuum  due  to  electronic 
excitation  in  the  metal.  This  excitation  (which  is  nonresonant 
with,  and  therefore  does  not  couple  to,  the  vibrational  modes)  is 
caused  by  the  visible  laser  used  in  SFG.  Antiresonance  would  still 
be  possible,  in  principle,  if  the  IR  laser  were  to  excite  an  electronic 
transition  in  the  SLG  at  the  frequency  of  a  C—H  stretching  mode. 
However,  no  coupling  was  observed,  which  is  consistent  with  the 
present  results.  In  the  experimental  configuration  used  in  ref  48, 
the  SLG  is  in  effect  a  metallic  layer  on  a  metallic  substrate,  and  the 
DDM  due  to  any  IR-active  electronic  excitation  within  the  SLG 
(or  the  Ir)  would  lie  parallel  to  the  surface45-47  and,  therefore, 
orthogonal  to  the  DMM  of  the  observed  C—H  modes.  It  is  also 
noted  that  the  observed  SFG  frequencies  (2563  and  2716  cm-1, 
respectively,  for  ortho  and  para  C—H  pairs)  are  much  lower  than 
the  frequencies  of  the  CH2  defect  sites  seen  here. 

The  effects  on  the  v(C— H)  spectrum  of  annealing  are  now 
considered  briefly.  Temperature-programmed  desorption  (TPD) 
studies49  of  H  on  carbon  nanosheets  show  several  H2  peaks  from 
~450  to  1 100  K,  indicating  a  variety  of  different  adsorption  sites 
and  bonding  configurations.  Because  similar  TPD  data  for  low 
coverages  of  H  on  SLG  are,  to  our  knowledge,  unavailable,  the 
nanosheet  results  are  used  as  a  guide.  Annealing  at  423  K  (i.e.,  just 
below  the  desorption  onset49)  in  flowing  Ar  had  no  significant 
effect  on  the  v(C— H)  spectrum  (not  shown).  However,  annealing 
at  673  K,  which  decomposes  some  nanosheet  sp3  C—H  sites, 
caused  major  changes  as  shown  in  Figure  6.  In  addition  to  the 
overall  loss  in  intensity,  indicating  removal  of  C—H  bonds,  there  is 
a  slight  shift  of  both  modes  to  higher  energy. 

A  major  effect  is  seen  in  vs(CH2),  which  now  exhibits  a 
strongly  antiresonant  (inverted)  line  shape  resembling  that  of 
vas(CH2).  This,  and  the  fact  that  vs(CH2)  is  independent  of  IR 
polarization  after  annealing  (not  shown),  indicates  that  the 
remaining  CH2  has  the  DDMs  of  both  v{C— H)  modes  lying 
in  the  SLG  plane.  This  accounts  for  the  lack  of  polarization 
dependence,  because  both  modes  cannowbe  excited  by  both  Epy 
and  Es>x  (Supporting  Information,  section  8),  and  for  the  strong 
coupling  of  both  modes  to  the  electronic  DDM.  Figure  6(inset) 
shows  the  implicit  structure  schematically.  An  in-plane  sp3 
(tetrahedral)  CH2  group  is  difficult  to  envision  in  SLG;  hence, 
we  are  led  to  consider  a  structure  like  a  terminal  C=CH2  group, 
which  might  occur  at  a  vacancy  or  domain  boundary  in  the  SLG 
or  at  the  periphery  of  the  small  BLG  patches  that  are  known  to 
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Figure  6.  The  effect  of  annealing  on  the  v(C— H)  spectrum.  The  data 
show  dR/R  (p-polarized)  for  CC14  exposure  at  P/P0  =  1.0  for  an  as- 
formed  sample  (a)  and  for  the  same  sample  after  a  673  K  vacuum  anneal 
(b).  Trace  (c)  shows  the  difference  (a  —  b),  and  spectra  have  displaced 
vertically  for  clarity.  The  inset  shows  schematically  a  CH2  configuration, 
viewed  along  the  z-axis  (i.e.,  the  surface  normal)  that  is  consistent  with 
spectrum  (b).  The  CH2  lies  mainly  in  the  SLG  plane  (gray),  as  do  the  vs 
(green)  and  vss  (red)  DDMs.  The  CH2  rotation  about  the  surface 
normal  is  arbitrary. 

exist  in  these  samples.  However,  in  organic  molecules,  v(C  —  H) 
for  such  species  typically  falls  above  ~3050  cm-1.  Hence,  the 
structure  after  annealing  cannot  yet  be  assigned  conclusively. 
There  may  also  be  an  effect  resulting  from  the  increased  contact9 
between  graphene  and  Si02  after  annealing. 

Again  using  the  nanosheet  H2  TPD  results49  as  a  guide,  we 
assume  that  the  species  remaining  after  the  673  K  anneal  are  also 
present,  in  the  same  chemical  state  and  in  the  same  concentration, 
before  annealing.  We  then  subtract  the  spectrum  after  annealing 
from  that  recorded  before  annealing  to  obtain  the  difference,  shown 
in  Figure  6c,  which  represents  the  spectrum  of  the  sp3  CH2  species 
removed  by  annealing.  This  is  qualitatively  similar  to  the  fitted  line 
shape  shown  in  Figure  5,  which  exhibits  a  “normal”  vs  and  a 
strongly  antiresonant  v ^  mode.  In  this  context,  small  sample-to- 
sample  variations  in  the  exact  C—H  line  shapes,  which  are  evident 
in  comparing  Figures  4—6,  can  be  understood  in  terms  of  varying 
relative  contributions  from  different  C—H  bonding  configurations. 

As  a  final  comment,  we  have  considered  (Supporting  Informa¬ 
tion,  section  12)  the  possibility  of  an  electrostatic  contribution  to 
the  observed  antiresonance  line  shapes  arising  from  a  coupling 
between  the  v(C— H)  dynamic  dipole  and  the  dipole  induced  in 
a  polarizable  adsorbate  such  as  CC14.  The  conclusion  is  that  such 
an  effect  cannot  account  for  the  observed  line  shapes. 

3.5.  Additional  Experiments.  In  an  effort  to  understand 
better  the  effects  reported  above,  we  now  discuss  three  further 
experiments  involving  (i)  analyte  desorption,  (ii)  multilayer  (>1 
layer)  graphene  films,  and  (iii)  the  effect  of  doping. 

Figure  7  shows  MIRS  data  versus  elapsed  time  during  N2 
purging  after  exposure  to  acetone  vapor  at  P/P0  =  1 .0.  Acetone  was 
selected  because  it  gives  an  intense  electronic  absorption  band  and 
because  the  v(C=0)  mode  is  sensitive  to  the  molecular 
environment.50  The  acetone  modes,  the  v(C— H)  features,  and 
a  large  part  of  the  electronic  absorption  band  all  vanish  within  A  f  = 
1  min,  which  is  the  maximum  attainable  time  resolution  in  the 
present  experiment.  The  remaining  features  are  a  weakened 
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Figure  7.  dR/R  (unpolarized)  versus  elapsed  time  (At)  during  purging 
with  pure  N2  after  exposure  to  acetone  vapor  at  P/P0  =  1.0.  Only  100 
scans  in  ~50  s  were  averaged,  and  At  is  the  elapsed  time  at  the  end  of  the 
data  collection.  The  At  =  0  spectrum  was  recorded  prior  to  shutting  off 
the  acetone  flow.  The  inset  shows  v(C=0)  at  two  different  times.  The 
feature  appearing  near  1450  cm  ,  with  nearly  constant  intensity  in  all 
spectra,  is  a  miscanceled  Si  multiphonon  absorption. 


electronic  absorption  band,  which  gradually  shifts  to  higher  energy 
with  continued  purging,  and  a  weak  v(C— O)  peak  that  slowly 
loses  intensity.  The  v(C— O)  feature  is  broad  and  red-shifted  to 
1700  cm-1  from  the  vapor-phase  value  of  1738  cm*1,  indicating50 
a  moderately  strong  interaction  with  a  range  of  polar  surface  sites, 
which  may  include  PMMA  residue  as  well  as  regions  of  bare  Si02 
(Supporting  Information,  section  6).  Even  for  At  as  short  as  1  min, 
there  is  no  sign  of  remnant  liquid  acetone  in  the  form  of  a  sharp 
l'(C=0)  feature  at  1715  cm*1.  Corresponding  data  obtained  as 
P/P0  is  increased  from  zero  (not  shown)  are  consistent  with  these 
results.  For  P/P0  below  '•“0.3,  the  weak  red-shifted  V’(C=0)  and 
the  weak  electronic  absorption  band  are  seen.  With  increasing 
P/ P0,  the  electronic  band  gains  intensity  and  shifts  to  lower  energy, 
while  the  acetone  modes  become  evident.  The  f(C— H)  anti¬ 
resonance  features  can  be  seen,  with  low  intensity,  at  P/ P0  as  low  as 
0.05  and  saturate  in  intensity  at  P/ P0  &  0.6. 

It  is  inferred  from  these  data  that  the  electronic  absorption 
results  from  an  essentially  local  effect.  At  high  P/Po,  the  whole 
SLG  layer  is  wetted,  including  the  nonpolar  C— H  sites,  which  do 
not  bond  strongly  to  acetone.  In  this  case,  the  electronic 
perturbation  extends  over  the  entire  SLG  layer,  and  the  anti¬ 
resonance  is  detected  clearly.  At  lower  P/Pq,  acetone  is  bound 
mainly  at  polar  sites,  and  the  electronic  perturbation  is  confined 
to  the  vicinity  of  these  sites,  which  probably  consist  of  a  low 
concentration  of  O-containing  defects.  The  C—H  sites  are,  at  this 
point,  uncovered  and  not  strongly  affected  by  acetone. 

The  slow  desorption  evident  in  Figure  7  might  also  indicate 
that,  at  high  P/Pq,  some  acetone  enters  the  space  between  SLG 
and  Si02  from  which  it  is  slow  to  escape  during  purging.  This  gap 
has  been  measured51  to  be  '-“8.2  A  for  as-prepared  SLG/ Si02  and 
is  assumed  to  be  filled  with  ambient  species  such  as  adsorbed 
H20.  However,  data  given  in  the  Supporting  Information, 
section  13,  show  that  acetone  does  not  enter  this  space,  although 
D20  vapor  easily  exchanges  with  the  H20  underlayer.  In  any 
case,  the  local  nature  of  the  effect  of  acetone  on  SLG  is  consistent 
with  the  proposal  that  it  is  related  to  adsorption-induced  strain. 


Figure  8.  Data  for  graphene  layers  of  varying  thickness,  (a)  SLG,  typical 
of  the  samples  used  in  the  present  work,  with  <5%  of  the  area  covered 
with  bilayer  regions;  (b)  an  SLG  layer  with  a  ~20%  of  the  area  covered 
with  bilayer  regions;  and  (c)  an  SLG  layer  with  '“40%  bilayer  regions.  All 
data  were  obtained  in  p-polarization  in  CC14  vapor  at  P/P0  =  1.0.  The 
dashed  line  shows  dR/R  =  0.  The  spectra  are  shown  as  recorded,  with  no 
displacement.  Spectrum  (a)  is  the  same  as  that  shown  in  Figure  2.  The 
weak,  unlabeled  feature  at  the  extreme  low-energy  end  of  (a)  and  (b)  is  a 
miscanceled  Si  multiphonon  absorption. 

The  above  discussion  suggests  a  reason  for  why  the  r(C— H) 
structure  (Figures  2  and  4)  is  relatively  weak  for  H20,  even 
though  an  electronic  absorption  is  clearly  evident.  If  regions  with 
sp5-hybridized  C—H  bonds  are  highly  hydrophobic,  they  would 
not  be  wetted  to  the  same  extent  as  the  rest  of  the  SLG,  even  at 
P/Pq  =  1.0,  and  would  remain  essentially  unaffected  by  the  H20. 
Pristine  regions  are  weakly  hydrophilic,  as  shown  by  the  disper¬ 
sion-corrected  DFT  results  described  above,  and  would  be 
wetted  at  elevated  P/Pq. 

Figure  8  shows  the  effects  of  increasing  the  graphene  thick¬ 
ness,  which  permits  an  assessment  of  the  importance  of  the  small 
BLG  regions  found  in  the  present  SLG  samples  (Supporting 
Information,  section  3)  and  also  facilitates  comparison  with 
previous  results.  The  fine  structure,  which  is  not  clear  in  Figure  8, 
is  shown  in  detail  in  the  Supporting  Information,  section  14. 
Several  points  are  noteworthy.  First,  the  C=C  stretching  mode 
at  about  1580  cm  1  is  IR-forbidden21  in  SLG  and  not  detectable 
in  the  present  SLG  data.  However,  it  is  IR-allowed  in  BLG  and 
FLG  and  is  clearly  observed  here  for  samples  with  increased 
bilayer  coverage.  This  feature  shows  a  clear  antiresonant  line 
shape,  as  in  previous  work, 24-27  although  the  sign  of  the 
asymmetry  differs  from  that  in  some  studies.24,27  The  f(C=0) 
mode  also  shows  a  clear  antiresonance.  In  both  cases,  the 
dynamic  dipole  lies  in  the  graphene  plane,  which  permits 
coupling  to  electronic  excitations  as  discussed  above.  Second, 
the  C—H  stretching  modes  become  weaker  with  increasing 
thickness.  This  is  consistent  with  previous  results24  in  which 
no  graphene  C—H  modes  were  seen  for  3—4-layer  films.  Because 
the  line  shapes  do  not  change  significantly,  one  concludes  that 
the  degree  of  coupling  to  the  electronic  continuum  is  essentially 
constant  and  that,  therefore,  the  loss  of  intensity  is  due  to  a  real 
decrease  in  the  number  of  C—H  bonds  per  unit  area. 

Third,  the  overall  shapes  of  the  broad-band  absorptions  for 
samples  with  significant  BLG  coverage  are  distinctly  different 
from  that  seen  for  SLG.  For  ~40%  BLG,  an  intense  and 
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Figure  9.  Data  showing  the  effect  of  a  strong,  hut  noncondensahle, 
dopant  (NH3)  on  a  ~20%  bilayer  sample.  Trace  (a)  shows,  for 
reference,  data  for  CC14  at  P/P0  =  1.0  as  in  previous  figures.  Traces 
(b)  — (d)  show  data  recorded,  after  completely  purging  the  CC14  with 
pure  N2,  in  a  100  mL/min  flow  of  0.1%  NH3  in  N2.  The  approximate 
total  exposure  times  are  (h)  4,  (c)  60,  and  (d)  220  min.  Trace  (e)  shows 
the  effect  of  purging  the  NH3:N2,  following  (d),  with  CCl4:N2atP/P0  =  1.0. 
The  spectra  were  recorded  in  p-polarization  and  have  not  been 
displaced  vertically.  The  dashed  line  shows  dR/R  =  0,  and  the  encircled, 
upward-pointing  hand  is  due  to  the  displacement  of  HzO  by  NH3  (see 
text). 

well-defined  band  is  seen  at  ~3030  cm-1  (0.38  eV)  with  a 
FWHM  of  "-'0.07  eV,  which  corresponds  closely  to  the  promi¬ 
nent  mid-IR  absorption  band  found23  in  transmission  spectra  of 
well-formed  BLG.  Simulation  (Supporting  Information,  section 
5)  shows  that  the  appearance  of  this  feature  in  dR/Ris  an  optical 
effect  resulting  from  the  influence  of  the  CC14  layer  on  the 
MIRS  “effective  thickness”  of  the  graphene  layer.  This  feature  is 
weak  in  the  ~20%  bilayer  sample,  for  which  the  second  layer 
is  more  fragmented  and  incomplete,  and  is  absent  in  SLG.  Thus, 
it  is  concluded  that  the  small  coverage  of  BLG  patches  found 
in  the  SLG  samples  has  no  significant  effect  on  the  present 
results. 

Finally,  the  data  in  Figure  9  show  that  doping  is  not,  in  itself,  a 
major  contributor  to  the  effects  observed  here.  For  another 
~20%  BLG  sample,  NH3  (an  n-type  dopant24)  leads  to  a  very 
intense  absorption  band  with  a  maximum  at  ''-'0.48—0.58  eV. 
After  a  total  exposure  of  ~220  min  to  0.1%  NFI3  in  N2,  the 
growth  of  this  band  had  slowed  but  not  yet  stopped.  The  ~20% 
BLG  film  used  here  is  similar  to  that  employed  in  Figure  8,  and 
the  CC14  spectrum  (Figure  9a)  is  also  similar  to  that  in  Figure  8a. 
The  antiresonant  v(C=C)  mode  gains  intensity  with  increasing 
NH3  coverage,24  and  a  broad,  inverted  band  is  seen  at 
~3430  cm  ,  which  is  ascribed  to  displacement  of  H20  from 
the  graphene/Si02  interface  by  NH3  (Supporting  Information, 
section  13).  The  NFI3  itself  is  not  detected,  presumably  due  to  a 
low  coverage  under  these  conditions.52  The  broad-band  feature 
seen  for  NFI3  is  different  from  that  observed  for  CC14  and  leads  to 
no  detectable  C— H  features.  Purging  the  NH3:N2  with  CC14:N2 
at  P/P0  =1.0  leads  to  the  immediate  appearance  of  the  C— H 
features  but  with  only  a  small  reduction  in  the  broad-band 
intensity.  Further  purging  with  saturated  CC14:N2  for  a  total  of 
45  min  (not  shown)  caused  only  a  small  additional  decrease  in 
the  broad-band  intensity. 


NIT3  bonds  only  very  weakly52  to  HOPG  and  does  not  form  a 
liquid-like  layer  at  room  temperature  under  the  present  condi¬ 
tions.  The  continued  growth  of  the  broad-band  absorption  over  a 
period  of  several  hours  in  a  constant  background  of  0.1%  NFI3: 
N2,  the  slow  reversal  of  the  NH3  effect  during  purging  (vs  the 
behavior  of  acetone,  Figure  7),  and  the  displacement  of  inter¬ 
facial  H20  all  indicate  that  NH3  diffuses  into  the  space  between 
graphene  and  Si02.  This  has  already  been  shown53  on  the  basis  of 
the  NH3  response  of  a  graphene/Si02  field-effect  transistor. 
Interfacial  NH3  is  able  to  bond  strongly54  to  Si—  OF!  via  an  O-H 
hydrogen  bond,  leaving  the  nonbonding  N  lone-pair  orbital  free 
to  donate  charge  to  the  graphene.  It  does  not,  however,  form  a 
“wetting”  layer  on  the  exposed  graphene  surface,  as  do  the  other 
small  molecules  studied  here,  and  does  not  affect  the  graphene 
electronic  structure  in  the  same  manner  as  do  these  other 
adsorbates. 

4.  CONCLUSIONS 

The  interaction  of  small  molecules  (CC14,  CS2,  H20,  and 
acetone)  with  SLG  has  been  observed  under  steady-state  condi¬ 
tions  using  IR  MIRS.  The  major  results,  most  of  which  have  not 
previously  been  reported,  are  as  follows. 

(1)  Abroad  and  intense  absorption  band,  spanning  the  range 
of  roughly  200—500  meV,  is  seen  as  a  result  of  exposure 
and  is  ascribed  to  an  electronic  excitation  resulting  from  a 
modification  of  the  SLG  band  structure  near  £P.  The 
qualitative  appearance  of  this  feature  and  its  dependence 
on  IR  polarization  vary  somewhat  from  sample  to 
sample,  which  suggests  an  extrinsic  effect  such  as  adsorp¬ 
tion-induced  strain.  However,  it  is  seen  for  all  samples 
studied  here  and  increases  in  intensity  as  CC14  <  H20  < 
acetone. 

(2)  In  an  effort  to  identify  the  source  of  the  electronic 
absorption  band,  ab  initio  DFT  calculations  have  been 
performed  to  model  the  interaction  of  H20  with  SLG 
with  and  without  an  Si02  substrate.  In  agreement  with 
previous  work,  proper  treatment  of  dispersion  (i.e.,  van 
der  Waals  interaction)  is  found  to  be  important  in 
obtaining  the  adsorption  energy.  An  adhesion  energy  of 
~1.5  eV/nm2  is  found  for  the  SLG/Si02  interface,  and 
the  configuration  of  the  adsorbed  H20  is  found  to  be 
vertical  for  free-standing  SLG  but  nearly  flat  for  SLG/ 
Si02.  It  is  found,  however,  that  the  electronic  absorption 
effects  seen  here  cannot  be  understood  in  terms  of  an 
ideal  H20/SLG/Si02  structure. 

(3)  Nonadiabatic  coupling  is  observed  between  C— H 
stretching  vibrations  and  the  absorption-induced  electro¬ 
nic  excitation.  This  leads  to  a  complex  antiresonance  line 
shape  for  the  vas(CH2),  and  to  some  extent  the  vs(CH2), 
mode  of  CH2  groups  present  as  defects  in  the  SLG.  On 
the  other  hand,  vibrational  modes  of  weakly  adsorbed 
molecular  species  do  not  undergo  this  coupling  and  thus 
exhibit  no  line  shape  anomalies.  Hence,  the  vibrating 
species  must  be  an  integral  part  of  the  SLG  lattice  for 
strong  coupling  to  occur.  A  linkage  exists  between  the 
strength  of  the  coupling  and  the  dependence  of  the  mode 
intensity  on  IR  polarization  because  both  are  related  to 
the  orientation  of  the  dynamic  dipole  moment  with 
respect  to  the  SLG  plane.  Thus,  the  effect  of  annealing 
on  the  polarization  and  line  shapes  of  the  v(C— H)  modes 
provides  indirect  evidence  of  restructuring  of  the  SLG. 
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These  results,  although  not  yet  fully  understood,  are  poten¬ 
tially  useful  in  several  respects.  The  effects  of  weakly  adsorbed 
small  molecules  on  the  SLG  electronic  structure  could  provide  a 
protocol  for  the  chemical  modification  of  SLG  electronic  proper¬ 
ties,  for  understanding  environmental  effects  on  these  properties, 
and  for  transduction  in  SLG  chemical  sensors.  The  vibrational 
modes  of  adsorbed  H,  as  seen  in  IR  spectroscopy,  have  long  been 
used  to  probe  semiconductor  surface  structure,  particularly 
under  steady-state  conditions  during  chemical  processing.  With 
the  added  dimension  of  nonadiabatic  coupling  to  the  adsorption- 
induced  electronic  excitation,  v(C  — H)  modes  are  potentially  of 
value  as  sensitive  indicators  of  changes  in  SLG  physical  and 
electronic  structure  during  treatment. 
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